arXiv:  1804.1 0395 v2  [astro-ph.CO]  18Jul2018 


Lethal  radiation  from  nearby  supernovae  helps  to  explain  the 

small  cosmological  constant 


Tomonori  Totani1,2,  Hidetoshi  Omiya1,  Takahiro  Sudoh1, 

Masakazu  A.R.  KobayashP,  Masahiro  Nagashima4 

department  of  Astronomy,  School  of  Science,  The  University  of  Tokyo, 

7-3-1  Hongo,  Bunkyo-ku,  Tokyo  113-0033,  Japan 
2Research  Center  for  the  Early  Universe,  School  of  Science,  The  University  of  Tokyo, 
7-3-1  Hongo,  Bunkyo-ku,  Tokyo  113-0033,  Japan 
3Faculty  of  Natural  Sciences,  National  Institute  of  Technology,  Kure  College, 
2-2-11  Agaminami,  Kure,  Hiroshima  737-8506,  Japan 
4Faculty  of  Education,  Bunkyo  University, 

3337  Minami-Ogishima,  Koshigaya-shi,  Saitama  343-8511,  Japan 


1 


Abstract 


The  observed  value  A0bs  of  the  cosmological  constant  A  is  extremely  smaller  than  theoretical 
expectations,  and  the  anthropic  argument  has  been  proposed  as  a  solution  to  this  problem  be¬ 
cause  galaxies  do  not  form  when  A  A0bs.  However,  the  contemporary  galaxy  formation 
theory  predicts  that  stars  form  even  with  a  high  value  of  A/A0bs  ~  50,  which  makes  the  an¬ 
thropic  argument  less  persuasive.  Here  we  calculate  the  probability  distribution  of  A  using  a 
model  of  cosmological  galaxy  formation,  considering  extinction  of  observers  caused  by  radi¬ 
ation  from  nearby  supernovae.  The  life  survival  probability  decreases  in  a  large  A  universe 
because  of  higher  stellar  density.  Using  a  reasonable  rate  of  lethal  supernovae,  we  find  that 
the  mean  expectation  value  of  A  can  be  close  to  Aobs,  and  hence  this  effect  may  be  essential  to 
understand  the  small  but  nonzero  value  of  A.  It  is  predicted  that  we  are  located  on  the  edge  of 
habitable  regions  about  stellar  density  in  the  Galaxy,  which  may  be  tested  by  future  exoplanet 
studies. 


1  Introduction 

The  concordance  cosmological  model  including  a  nonzero  cosmological  constant,  A,  is  ac¬ 
cepted  as  the  best  description  of  the  universe  observed  (Weinberg  et  al.  2013;  Bull  et  al.  2016). 
However,  the  observed  value  A0bs  is  smaller  than  theoretical  expectations  by  at  least  50  orders 
of  magnitude.  Furthermore,  A  works  as  the  constant  energy  density  of  vacuum,  and  its  value  is 
curiously  similar  to  the  decreasing  matter  energy  density  just  now  in  the  long  history  of  the  uni¬ 
verse.  The  extreme  smallness  and  coincidence  are  perhaps  the  most  difficult  problem  in  physics 
(Carroll  2001;  Caldwell  and  Kamionkowski  2009).  Lack  of  persuasive  explanations  has  moti¬ 
vated  to  consider  the  anthropic  argument,  assuming  that  A  is  determined  by  a  stochastic  process 
in  the  early  universe.  An  observer  will  not  appear  in  a  universe  of  A  Aobs,  because  structure 
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formation  is  suppressed  before  galaxies  form  once  the  cosmic  expansion  is  accelerated  by  A 


(Barrow  and  Tipler  1986;  Weinberg  1987;  Efstathiou  1995;  Martel  et  al.  1998;  Garriga  et  al. 
2000;  Peacock  200711. 

A  wide  range  of  theoretical  possibilities  about  stochastic  determination  of  A  have  been 
discussed  (e.g.  Polchinski  2006;  Totani  2017).  Here  we  consider  the  simplest  version  in  which 
A  is  the  only  random  variable  while  other  physical  constants  and  the  total  amount  of  matter 
in  the  universe  are  unchanged.  Though  some  theories  predict  that  not  only  A  but  also  other 
physical  constants  may  change,  we  leave  these  possibilities  for  future  studies.  If  the  prior 
probability  distribution  dPp/dA  of  A  extends  down  to  A  =  0  with  a  nonzero  value,  a  natural 
distribution  is  flat  per  unit  A  (i.e.  dPp/dA  =  const.  1]  around  A  =  0,  because  physically  natural 
scales  of  A  are  much  larger  than  Aobs.  In  this  work  we  assume  this  prior  distribution  following 
previous  studies  (e.g.,  Efstathiou  1995;  Martel  et  al.  1998;  Garriga  et  al.  2000;  Peacock  2007). 
Then  the  probability  distribution  of  A  for  an  observer,  dPa/dA,  should  be  proportional  to  the 
number  of  observers  per  unit  mass  of  cosmic  matter,  n( A).  A  very  small  value  of  |A|  <C  Aobs 
is  statistically  unlikely,  and  hence  the  coincidence  is  also  explained.  Though  A  can  also  be 
negative,  the  cosmic  expansion  would  be  pulled  back  to  collapse  before  an  age  of  ~10  Gyr 
when  A  £  —  Aobs,  allowing  no  observer  like  us  in  a  universe  with  A  -C  —  Aobs.  This  means  that 
inclusion  of  negative  A  does  not  significantly  affect  the  anthropic  argument,  and  we  consider 
the  A  distribution  only  in  A  >  0. 

More  quantitatively,  n( A)  can  be  calculated  based  on  the  theory  of  galaxy  formation,  if  we 
assume  that  the  number  of  observers  is  proportional  to  stellar  mass  produced  in  the  universe 
(Sudoh  et  al.  2017;  Barnes  et  al.  2018).  The  calculation  of  Sudoh  et  al.  (2017)  by  a  semi- 


1  See  also  Lineweaver  &  Egan  (2007)  and  Egan  and  Lineweaver  (2008)  for  anthropic  expectation  about  the  ratio 
of  matter  to  vacuum  density  of  the  universe,  under  the  condition  that  the  cosmological  constant  is  the  observed 
value  (i.e.,  A  =  Aobs). 

2If  A  is  bound  to  be  positive,  a  flat  distribution  per  unit  In  A  may  also  be  possible  without  any  characteristic 
scale.  However,  in  this  case  we  need  to  introduce  a  lower  cut-off  at  A  <C  Aobs,  otherwise  the  distribution  is 
mathematically  ill-defined.  We  do  not  know  any  low-energy  physics  leading  to  such  a  cut-off. 
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analytic  model  of  cosmological  galaxy  formation  (Nagashima  and  Yoshii  2004)  is  shown  by 
solid  curves  in  Fig.  [H  where  n( A)  is  calculated  by  stellar  masqj  produced  up  to  a  fixed  age  of 
the  universe  (15  Gyr),  and  the  cosmological  parameters  except  for  A  are  fixed  to  the  2015  deter¬ 
mination  by  Planck  (Planck  Collaboration  2016).  This  model  is  based  on  the  standard  picture 
of  hierarchical  structure  formation  driven  by  cold  dark  matter.  Dark  matter  halos  are  generated 
by  the  Monte-Carlo  method,  so  that  their  formation  rate  and  merger  history  are  consistent  with 
the  structure  formation  theory.  Then  galaxies  grow  in  dark  halos  by  star  formation,  which  are 
calculated  considering  baryonic  physics  such  as  gas  cooling,  supernova  feedback,  and  galaxy 
mergers.  Various  observations,  including  galaxy  luminosity  functions  of  local  as  well  as  high- 
redshift  galaxies,  are  broadly  consistent  with  the  predictions  by  this  model  (Nagashima  and 
Yoshii  2004;  Kashikawa  et  al.  2006;  Kobayashi  et  al.  2007).  The  predicted  distribution  of 
dPQ/d{ In  A)  peaks  at  A  ~  20  A0bs,  and  the  distribution  extends  to  A  ~  100  A0bs.  The  proba¬ 
bility  of  finding  A  <  Aobs  is  only  PQ( A  <  Aobs)  =  6.7%.  More  recently,  Barnes  et  al.  (2018) 
found  the  distribution  peak  at  A/Aobs  =50-60  and  Pa{ A  <  Aobs)  =  2%  using  a  numerical  sim¬ 
ulation  of  galaxy  formation,  which  is  roughly  consistent  with  Sudoh  et  al.  considering  model 
uncertainties.  This  may  imply  that  the  value  of  A  is  not  determined  by  the  anthropic  mechanism. 

However,  the  amount  of  star  formation  is  not  the  only  factor  to  determine  the  number  of 
observers,  but  we  should  also  consider  planet  habitability  around  a  star,  which  depends  on  loca¬ 
tions  and  environments  in  a  galaxy  (Gonzalez  et  al.  2001;  Lineweaver  et  al.  2004;  Gowanlock 
et  al.  2011;  Forgan  et  al.  2017).  Radiation  from  nearby  supemovae  is  widely  discussed  as  a 
possible  cause  of  extinction  of  terrestrial  organisms  on  a  planet,  because  high  energy  radiation 
(gamma-rays  and  cosmic-rays)  from  a  core  collapse  supernova  within  ~10  pc  would  have  sig¬ 
nificant  effects  on  the  ozone  layer  (Ruderman  1974;  Whitten  et  al.  1976;  Reid  et  al.  1978; 


3Sudoh  et  al.  (2017)  also  considered  a  dependence  of  n( A)  on  metallicity,  but  found  that  this  effect  is  actually 
small  (see  their  Fig.  4).  This  is  because  low  metallicity  galaxies  are  generally  dwarf  galaxies,  and  star  formation  in 
such  galaxies  is  suppressed  by  supernova  feedback,  making  their  contribution  to  cosmic  stellar  mass  budget  small. 
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Gehrels  et  al.  2003).  Interestingly,  the  expected  number  of  such  nearby  supernovae  in  the  so¬ 
lar  neighborhood  is  of  order  unity  during  the  time  span  of  terrestrial  organisms  (~  0.5  Gyr)  in 
the  history  of  Earth  (Gehrels  et  al.  2003).  This  coincidence  implies  that  we  may  be  living  on 
the  edge  of  habitable  regions  concerning  stellar  density,  and  an  intelligent  observer  may  not  be 
allowed  to  appear  in  higher  density  regions.  In  this  work  we  make  a  quantitative  assessment  of 
this  effect  on  the  probability  distribution  of  A. 

2  Stellar  Densities  in  Galaxy  Formation 

The  key  quantity  in  galaxy  formation  to  determine  stellar  density  is  the  internal  density  of  a 
gravitationally  collapsed  dark  halo.  The  physics  of  dark  halo  formation  can  be  understood  by  the 
analytic  spherical  collapse  model  (see,  e.g.  Cooray  and  Sheth  2002  for  a  review).  The  internal 
density  (virial  density)  after  the  collapse,  pv ir,  does  not  depend  on  halo  masses  but  it  decreases 
with  cosmic  time,  as  shown  in  Fig.  [2]for  various  values  of  A.  When  A  is  large,  decreasing  energy 
density  of  matter  becomes  lower  than  A  in  an  earlier  epoch.  After  this  transition,  pv]T  becomes 
constant,  and  the  number  of  collapsing  halos  rapidly  decreases.  Therefore,  in  a  universe  with 
A  =  50  Aobs,  internal  density  of  any  dark  halo  is  more  than  10  times  higher  than  that  of  a  halo 
collapsing  at  a  cosmic  time  of  10  Gyr  in  the  universe  that  we  observe.  Stars  are  expected  to  form 
only  after  the  hot  diffuse  gas  in  a  halo  cools  and  contracts,  but  the  final  size  of  a  rotationally 
supported  gas  disk  is  proportional  to  the  virial  radius  of  the  host  dark  halo,  because  the  mean 
specific  angular  momentum  given  to  halo  gas  is  approximately  universal  (Mo  et  al.  1998). 
Therefore  pvir  is  a  good  indicator  of  stellar  density  p*  in  a  halo. 

The  Sun  formed  9.2  Gyr  after  the  big  bang  (a  redshift  of  z  =  0.42).  According  to  the  cosmic 
star  formation  history,  about  90%  of  all  stars  formed  by  the  present  time  (13.8  Gyr)  are  older 
than  the  Sun  (Borch  et  al.  2006),  and  hence  in  higher  density  regions.  This  also  implies  the 
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effect  of  supernova  radiation  on  life;  if  this  effect  is  not  working,  it  is  statistically  more  likely 
to  find  ourselves  in  higher  density  regions  than  the  solar  neighborhood.  It  should  also  be  noted 
that  such  regions  have  more  time  for  evolution  of  life. 

For  the  first  simple  calculation  of  n(A)  with  the  nearby  supernova  effect,  we  estimate  stellar 
densities  in  the  galaxy  formation  model  as  follows.  Any  model  galaxy  is  either  the  central 
galaxy  in  a  dark  halo  or  a  satellite  galaxy  trapped  by  a  dark  halo  as  a  result  of  mergers  with 
smaller  halos.  The  density  of  a  stellar  population  born  in  a  central  galaxy  is  assumed  to  be 
p*  oc  pvir  of  the  host  halo,  while  that  in  a  satellite  galaxy  is  estimated  from  pvir  of  the  last  host 
halo  in  which  the  galaxy  existed  as  the  central  galaxy.  In  some  cases  stellar  density  would  be 
changed  by  a  galaxy-galaxy  merger  in  a  dark  halo,  but  we  ignore  this  effect  for  simplicity. 

The  halo  virial  density  of  the  Milky  Way  (MW)  Galaxy  is  estimated  as  pvir,MW  ~  1-4  x 
104  Me  kpc-3  (Klypin  et  al.  2002),  which  can  be  used  as  an  indicator  of  p*  around  the  solar 
neighborhood  (p*,0),  assuming  that  the  Sun  is  a  typical  star  formed  around  the  main  collapse 
epoch  of  the  present  MW  halo.  Then  we  can  determine  the  stellar  density  relative  to  the  solar 
neighborhood  value  as  p*/p*,©  =  pvir/Pvir,MW>  for  all  stellar  populations  formed  at  each  time 
step  in  the  galaxy  formation  model.  It  would  be  possible  to  introduce  more  detailed  modelings 
about  internal  structures  of  a  galaxy,  but  they  depend  on  complicated  baryon  physics  and  model 
uncertainties  would  become  large.  The  primary  cosmological  effect  of  A  to  make  stellar  density 
higher  has  reasonably  been  taken  into  account  by  this  treatment. 

3  Probability  Distribution  of  A  under  the  Nearby  Supernova 
Effect 

Let  iVexp  for  a  stellar  population  be  the  expected  number  of  sterilizing  events  during  the  time 
for  evolution  of  life  to  an  observer.  We  introduce  a  parameter  NexPt&,  the  value  around  the  solar 
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neighborhood,  which  is  rather  uncertain  but  of  order  unity  if  we  consider  nearby  supernovae. 
Obviously  Nexp  should  increase  with  stellar  density,  i.e.,  Nexp  =  (p*/P*,o)  AfexP)0.  Here,  we 
implicitly  assumed  that  the  supernova  rate  is  the  same  for  all  stellar  populations,  though  the  rate 
of  core-collapse  supernovae,  which  are  related  to  young  stellar  populations,  should  be  small  in 
old  stellar  populations.  However,  type  la  supemovae  (SNe  la)  would  also  affect  habitability, 
which  are  related  to  white  dwarfs  and  occur  also  in  old  stellar  populations  (Totani  et  al.  2008). 
The  SN  la  rate  is  about  ten  times  lower  than  the  core-collapse  rate  in  MW  Galaxy,  but  they 
may  have  more  significant  effect  on  habitability  because  of  higher  luminosity  (Gowanlock  et  al. 
2011).  Lethal  events  occurring  in  old  stellar  populations  are  important  in  galaxies  dominated 
by  such  populations,  like  bulges  or  elliptical  galaxies. 

From  the  Poisson  statistics,  the  probability  of  survival  with  no  lethal  event  is  Ps  =  exjp(—Nexp). 
Then  n(A)  is  calculated  by  integrating  all  stellar  mass  produced  in  the  universe,  but  with  the 
weight  of  Ps  for  each  stellar  population  in  a  galaxy.  Figure  Q]  shows  n(A)  and  the  probability 
distributions  of  A  for  some  different  values  of  NexPtQ.  The  location  of  the  distribution  peak  is 
reduced  to  A/Aobs  ~  4  when  NexPte  =  1,  and  P0(A  <  Aobs)  is  increased  to  19%.  If  A^xp,©  =  3, 
a  rather  small  probability  of  Ps(NexPjQ )  =5%  means  that  we  are  a  lucky  survivor,  but  the 
distribution  peak  is  further  reduced  to  A/Aobs  ~  2  with  the  probability  of  PC(A  <  Aobs)  = 
41%. 


4  Discussion 

There  are  still  many  uncertain  aspects  about  this  calculation.  The  critical  distance  to  a  lethal 
supernova  suffers  from  uncertainties  about  the  effect  of  high  energy  particles  on  the  Earth  at¬ 
mosphere,  and  the  degree  of  damage  by  ultraviolet  light  on  organisms.  Stellar  density  should 
be  calculated  by  simulations  resolving  internal  structure  of  galaxies  in  future  studies,  and  su- 
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pernova  rates  should  be  calculated  by  local  star  formation  history  in  a  galaxy.  However,  the 
calculation  shown  here  is  based  on  reasonable  parameters  and  treatments,  and  hence  lethal  radi¬ 
ation  from  nearby  supernovae  must  be  included  in  the  anthropic  consideration  about  A,  which 
may  be  the  only  way  proposed  so  far  to  understand  both  the  extreme  smallness  and  coincidence. 
Supemovae  played  an  important  role  for  humankind  to  discover  the  nonzero  cosmological  con¬ 
stant  (Riess  et  al.  1998;  Perlmutter  et  al.  1999),  but  they  may  also  be  responsible  for  its  small 
value. 

If  there  is  any  physical  effect  other  than  supernovae  to  reduce  habitability  in  high  stellar 
density  regions,  it  would  also  contribute  to  making  the  expected  value  of  Aobs  smaller  (e.g., 
comet  bombardment  by  a  field  star  passage  or  very  wide  binary  star  systems  affected  by  the 
Galactic  potential,  see  Kaib  2018).  Especially,  gamma-ray  bursts  (GRBs)  have  been  discussed 
as  a  hazardous  astronomical  event,  whose  impact  may  be  comparable  or  even  larger  than  super¬ 
novae  (Thorsett  1995;  Melott  et  al.  2004;  Piran  and  Jimenez  2014).  GRB  rate  is  much  lower 
than  supernova  rates,  but  their  luminosity  is  much  brighter,  and  hence  the  maximum  distance 
to  a  lethal  GRB  is  comparable  to  a  typical  size  of  a  galaxy.  Piran  et  al.  (2016)  proposed  that 
lethal  radiation  from  GRBs  in  nearby  dwarf  galaxies  prohibits  an  observer  to  exist  in  a  low  A 
universe  (0  <  A  <C  Aobs),  which  would  be  important  if  the  prior  A  distribution  is  weighted  to¬ 
wards  smaller  values,  rather  than  being  uniform  per  unit  A.  It  should  be  noted  that  long-duration 
GRBs  (LGRBs)  are  related  only  to  young  stellar  populations,  but  lethal  events  in  old  stellar  pop¬ 
ulations  are  also  important  to  reduce  Aobs  by  the  effect  considered  in  our  work.  Short-duration 
GRBs  occur  also  in  old  stellar  populations,  like  SNe  la,  but  their  effect  on  habitability  is  likely 
much  smaller  than  LGRBs  due  to  the  smaller  energy  emitted  per  event.  Another  important  dif¬ 
ference  from  supernovae  is  that  LGRBs  occur  preferentially  in  low  metallicity  regions.  Since 
high  stellar  density  regions  are  generally  metal  rich  in  a  galaxy,  this  should  weaken  the  effect 
of  LGRBs  on  habitability  in  high  density  regions. 
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If  the  nearby  supernova  effect  is  indeed  responsible  for  the  observed  small  value  of  A  as 
proposed  here,  a  prediction  is  that  we  should  be  located  near  the  edge  of  habitable  regions  about 
stellar  density  in  a  galaxy.  This  may  be  tested  in  future  by  development  of  observational  studies 
on  exoplanets  and  their  habitability.  Exoplanets  in  regions  of  higher  stellar  density  would  show 
smaller  probability  of  biomarker  detection,  even  if  they  are  apparently  habitable. 
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Figure  1:  Top  panel:  the  number  of  observers  n( A)  per  unit  mass  of  cosmic  matter,  normalized 
to  unity  in  the  limit  of  A  =  0,  for  several  different  values  of  Nex P:0  (the  expected  number  of 
lethal  events  in  the  solar  neighborhood  during  the  time  for  evolution  of  life  to  an  observer). 
When  there  is  no  supernova  effect  (i.e.,  NexP:Q  =  0),  n(A)  is  proportional  to  stellar  mass  pro¬ 
duced  in  the  universe  up  to  a  cosmic  age  of  15  Gyr.  Bottom  panel:  the  probability  distribution 
of  A  for  an  observer  per  unit  In  A,  dPQ/d(  In  A)  =  AdPa/dA  oc  An(A),  assuming  a  flat  prior 
distribution  of  dPp/dA  =  const. 
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Figure  2:  Virial  density  Mvir/ (47rr;3ir/3)  of  a  dark  halo  as  a  function  of  its  collapsing  time  for 
different  values  of  A,  where  Mvir  and  rvir  are  the  virial  mass  (baryon  plus  dark  matter)  and  virial 
radius,  respectively. 
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